Critical properties of the iV-color London model 
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The critical properties of iV-color London model are studied in d = 2 + 1 dimensions. The model 
is dualized to a theory of N vortex fields interacting through a Coulomb and a screened potential. 
The model with N = 2 shows two anomalies in the specific heat. From the critical exponents a and 
v, the mass of the gauge field, and the vortex correlation functions, we conclude that one anomaly 
corresponds to an inverted 3Dxy fixed point, while the other corresponds to a 3Dxy fixed point. 
There are N fixed points, namely one corresponding to an inverted 3Dxy fixed point, and N — 1 
corresponding to neutral SDxy fixed points. This represents a novel type of quantum fluid, where 
superfluid modes arise out of charged condensates. 

PACS numbers: 11.15.Ha, 11.27.+d, 71.10.Hf, 74.10.+V, 74.90. +n 



Ginzburg-Landau (GL) theories with several complex 
scalar matter fields minimally coupled to one gauge field 
are of interest in a wide variety of systems, such as multi- 
ple component (color) superconductors, metallic phases 
of light atoms such as hydrogen 0, and as effec- 
tive theories for easy-plane quantum antiferromagnets 
[3- ■ The model also is highly relevant in particle 
physics where it is called two-Higgs doublet model . In 
metallic hydrogen the scalar fields represent Cooper pairs 
of electrons and protons, which excludes the possibility 
of inter-color pair tunneling, i.e. there is no Josephson 
coupling between different components of the condensate. 
The same two-color action in (2 + 1) dimensions, where 
the matter fields originate in a bosonic representation of 
spin operators, is claimed to be the critical sector of a 
field theory separating a Neel state and a paramagnetic 
(valence bond ordered) state of a two dimensional quan- 
tum antiferromagnet at zero temperature with an easy- 
plane anisotropy present . This happens because, 
although the effective description of the antiferromag- 
net involves an a priori compact gauge field, it must be 
supplemented by Berry-phase terms in order to properly 
describe S — 1/2 spin systems 0,0. Berry-phases cancel 
the effects of monopoles at the critical point 00. In this 
paper, we point out novel physics of the quantum fluid 
that arises out of an TV-color charged condensate when 
no intercolor Josephson coupling is present. 

For a detailed analysis of the phase transitions in such 
a generalized GL model, we study an N component GL 
theory in (2 + 1) dimensions with no Josephson coupling 
term. The model is defined by N complex scalar fields 
{\J/( Q )(r) | a — 1 . . . N} coupled through the charge e to 
a fluctuating gauge field A(r), with Hamiltonian 



H 
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where is the a-component condensate mass. The 

potential y({^ a )(r)}) is assumed to be only a function 
of ^(^(r)! 2 . The model is studied in the phase-only 



(London) approximation \lK a )(r) = | exp[i6 

and is discretized on a lattice with spacing a = 
In the Villain approximation the partition function reads 
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where n( Q )(r) are integer vector fields ensuring 27T pe- 
riodicity, and the lattice position index vector r of the 
fields is suppressed. The symbol A denotes the lattice 
difference operator and j3 = 1/T is the inverse tempera- 
ture. Here, we stress the importance of keeping track of 
the 2tt periodicity of the individual phases. The kinetic 
energy terms are linearized by introducing N auxiliary 
fields v( Q ). Integration over all 6^ produces the local 
constraints A-v^) = 0, which are fulfilled by the replace- 
ment v( Q ) — > A x h^ a K We recognize hW as the dual 
gauge fields of the theory. By fixing the gauge ni = 
and performing a partial integration we may introduce 
the vortex fields mS°^ — A x nW. We integrate out the 
gauge field A and get a theory in the dual gauge fields 
h( a ) and the vortex fields m( Q ) where A • m( Q ) = 
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where \^\ 2 = |^ a) | 2 /MW. Note how the algebraic 
sum of the dual photon fields is massive. This differs 
from the case N = 1, where e produces one massive dual 
photon with bare mass e 2 /2, and the model describes a 



2 



vortex field m interacting through a massive dual gauge 
field h. However, when N > 2, since A • m'"^ = 0, a 
gauge transformation hW — * h^+Ag^ for a — 1 ... N 
leaves the action invariant if one of the gauge fields, say 
h^- 1 compensates the sum in the last term in (0 with 

Integrating out the dual gauge fields we get a general- 
ized theory of N interacting vortex fields 
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where c^j, is the Kronecker-delta 



and the vortex inter- 
action potential D^ a ' r, \r) is the inverse discrete Fourier 
transform of D^ '^ (q), where 
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Here, 

= e 2, 2 is the square of the bare inverse screen- 
ing length in the intervortex interaction, and |Q q | 2 is 
the Fourier representation of the lattice Laplace opera- 
tor. The first term of the vortex interaction potential 
is a Yukava screened potential, while the second term 
mediates long range Coulomb interaction between vortex 
fields. If N = 1 the latter cancels out exactly and we are 
left with the well studied vortex theory of the GL model 
which has a charged fixed point for e / [ToL [i"0| . For 
N > 2 wc find a theory of vortex loops of N colors in- 
teracting through long range Coulomb interaction. If N 
grows to infinity, ip 2 — > oo and the vortex fields interact 
via a diagonal unscreened N x N Coulomb matrix. This 
reflects the inability of one single gauge field A to screen 
a large number of vortex species. The case N > 2 has 
features with no counterpart in the case N = 1 "0, [7^/ . 
namely neutral superfluid modes arising out of charged 
condensates. 

The above vortex system may be formulated as a field 
theory, introducing N complex matter fields (f>^ for each 
vortex species, minimally coupled to the dual gauge fields 
h("). This generalizes the dual theory for N = 1 pio- 



neered in The theory reads (see also |5j|) 
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Here, we have added chemical potential (core-energy) 
terms for the vortices, as well as steric short-range repul- 
sion interactions between vortex elements. In the N = 1 



case, a RG treatment of the mass term of the dual gauge 
field yields de 2 /dlnl = e 2 , and hence this term scales 
up, suppressing the dual gauge field. Correspondingly, 
for N > 2, this suppresses ^2 a ^ a \ but not each indi- 
vidual dual gauge field. For the particular case N = 2, 
assuming the same to hold, we end up with a gauge the- 
ory of two complex matter fields coupled minimally to 
one massless gauge field, which was also precisely the 
starting point. Thus the theory is self-dual for N = 2 
0,0. For N = 1, it is known that a charged theory in 
d = 2 + 1 dualizes into a |0| 4 theory and vice versa \U\ . 
The vortex tangle of the SDxy model is incompressible 
and the dual theory is a massless gauge theory such that 
((f)) ^ is prohibited. For e ^ 0, the dual theory has 
global symmetry, and vortex condensation and (<f>) ^ 
is possible |ll|. 

For N — 2, Monte Carlo (MC) simulations have been 
carried out for the action (0J with parameters \tp^ | 2 = 



1/2, 



2)12 



1, e 2 = 1/4, and m 2 , = 3/8. Here, \ip { 
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and \ip( 2 ^ \ 2 have been chosen to have well-separated bare 
energy scales associated with the twist of the two types of 
phases, and mo has been chosen to be of the order of the 
inverse lattice spacing in the problem to avoid difficult 
finite-size effects. One MC update consists of inserting 
elementary vortex loops of random direction and species 
according to the Metropolis algorithm. 

We observe two anomalies in the specific heat at T c \ 
and T C 2 where T c \ < T C 2 . We find T c \ and T C 2 from scal- 
ing of the second moment of the action ((SV — {Sv)) 2 ) 
to be T c i = 1.4(6) and T c2 = 2.7(8). To check the criti- 
cality of these anomalies we have calculated the critical 
exponents a and v by applying finite size scaling (FSS) 
of M 3 = {(S v - (Sv)) 3 ) E3- The peak to peak value 
of this quantity scales with system size L as 
the width between the peaks scales as L^ 1 ^ . The ad- 
vantage of this is that asymptotically correct behavior 
is reached for practical system sizes. The FSS plots for 
system sizes L — 4,6, 8, 10, 12, 14, 16, 20, 24 are shown in 
Fig. ^ From the scaling we conclude that both anomalies 
are in fact critical points, and we obtain a = — 0.02±0.02 
and v = 0.67 ± 0.01 for T cl and a = -0.03 ± 0.02 and 
v = 0.67±0.01 for T C 2. These values are consistent with 
those of the SDxy and the inverted 3Dxy universality 
classes found with high precision to be a = —0.0146(8) 
and v = 0.67155(3) 

To characterize these phase transitions further, 
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The masses of A and E Q ^ are defined by 
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FIG. 1: The FSS of the peak to peak value of the third mo- 
ment AM3 labeled (□) and (+) for T c \ and T C 2 respectively. 
The scaling of the width between the peaks A/3 labeled (a) 
and ( x ) for T c \ and T C 2 respectively. The lines are power law 
fits to the data for L > 6 used to extract a and v. 



lim ? ^o 2g A (q)- 1 /(3 and m| h = lim^o 2/?V> 2 £sh(?)- 1 . 

We briefly review the case N = 1 [llj. The dual 
field theory of the neutral fixed point (toq = 0) is a 
charged theory describing an incompressible vortex tan- 
gle. The leading behavior of the vortex correlator is 



lim^ 27r 2 /3G( +) (<7) ~ [l-C 2 (T)]q 2 



C 3 (T)q 2 +^, 



and q 2 + C 4 (T)q A for T < T c , T = T c , and T > T c respec- 
tively. For T < T c we have m|. h = (N =1), however for 
T > T c the l/q 2 terms in Q^h (<?) cancel out exactly and 
this mass attains an expectation value. At the charged 
fixed point (to 2 , ^ 0) of the GL model, the effective field 
theory of the vortices is a neutral theory. The vortex 
tangle is compressible with a scaling ansatz for the vor- 
tex correlator lim 9 ^ G (+) (q) ~ q 2 , q 2 ~ VA , and c(T) for 
T < T Cl T = T c , and T > T c , respectively. Consequently, 
from J7J), the mass to a drops to zero at T c , and the mass 
of the dual gauge field mh is finite for all temperatures 
and has a kink at T c . Renormalization group ar gum ents 
yield 77 a = 4 — d where d is the dimensionality Jlj^- Il2l | , 
which has recently been verified numerically [111 |l6j . 
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FIG. 2: G (+) (q) for N = 2, L = 32. For T = 2.86 > T c2 , 
T = 2.76 ~ T c2 , and T = 2.63 < T c2 , lim q ^ G {+) {q) ~ c(T), 
~ and ~ g 2 , respectively. 



The vortex correlator for N = 2 is sampled in real 
space and G^ (q) is found by discrete Fourier transfor- 
mation, it is shown in Fig. (J2J. At T — T cl the leading 
behavior is G^ + \q) ~ g 2 on both sides of T c i. Conse- 
quently, due to 0, to a and m^h are finite in this regime. 
This shows that the vortex tangle is incompressible and 
that the anomalous scaling dimension t/a = 0, which 
corresponds to a neutral fixed point. Below T c2 the dom- 
inant behavior is G^ + '(q) ~ q 2 whereas G^(q) ~ c(T) 
above T c2 . At T = T c2 , G i+) (q) ~ q indicating t]a = 1. 
Accordingly, to a is finite below T C 2 and zero for T > T C 2. 

For T < T c2 , mA scales according to Qxiq)^ 1 ^ — 
m\ + Cq 2 ~ r,A + 0(q s ) for small q where 6 > 2 - r/A 
|lfl| . with a corresponding Ansatz for <?sh(<?)- For each 
coupling we fit 5a(q , ) _1 data from system sizes L = 
8, 12, 20, 32 to the Ansatz. The results for toa (and ma, 
found similarly), are given in Fig. |3J The system ex- 
hibits Higgs mechanism at T = T c2 when toa drops to 
zero. Furthermore to a has a kink at T c \ due to ordering 
of The anomalies in tua coincide precisely with T c \ 
and T c2 determined from scaling of {{Sy — (SV)) 2 )- Note 
also how msh changes abruptly at T c2 . This is due to a 
sudden change in screening of Y] 1 h^"- 1 by the vortex- 
loop proliferation at T = T c2 , giving an abrupt increase 
in msh, analogously to what happens for N = l,e ^ 
|ll| . Above T c2 , A is massless, giving a compressible vor- 




FIG. 3: The mass mA (♦) and 1 — mo/fish (+) found from 
Eq. 0. Two non-analyticities can be seen in tua at T c i and 
T C 2, corresponding a neutral fixed point and a charged Higgs 
fixed point, respectively. An abrupt increase in m^h due to 
vortex condensation is located at T C 2. 

tex tangle which accesses configurational entropy better 
than an incompressible one. Below T c2l A is massive and 
merely renormalizes |\&| 4 terms in Eq. QJ. The theory is 
effectively a l^l 4 theory in this regime. Thus, the remain- 
ing proliferated vortex species originating in the matter 
fields with lower bare stiffnesses form vortex tangles as if 
they originated in a neutral superfluid. For the general 
N Higgs mass is generated at the highest critical 

temperature, after which A merely renormalizes the \^>\ 4 
term, such that the Higgs fixed point is followed by N — 1 
neutral fixed points as the temperature is lowered. 
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We now discuss the vortex mode — m^ 2 \ demon- 
strating that it should be identified as a superfluid 
mode in the system. Its properties are controlled by 
<?Ah(<?) = (Ihq 1 ^ — h^| 2 }. A dual Higgs phenomenon for 
N = 2,T = T c \ involving t/Ah(?) may be demonstrated 
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The G^(q) correlation function is always ~ q 2 ,q — > 0, 
but has a nonanalytic coefficient of q 2 , determined by 
the helicity modulus T of the neutral mode — m^ 2 ). 
When T vanishes at T c \ through a disordering of 
thus destroying the superfluid neutral mode, the first 
and second term in the bracket cancel, which in turn 
cancels the 1/q 2 term in GAh(q)- This produces a dual 
Higgs mass mAh defined by 0Ah(?) ~ l/(<? 2 + m 2 Ah ) for 
T > T c i. The remaining terms in Eq. JHJ contribute 
to determining the actual value of mAh- Thus, while 
h (i) + h (2) ig always ma ssive, c f. Eq. ©, hM - h^ 2 ' 
is massless below T c \ and massive above T c \. Therefore 
_ h( 2 ) plays the role of a gauge degree of freedom, 
providing a dual counterpart to A in Eq. Q. This is 
evident when |i/> (1) | 2 = |i/> (2) | 2 . Then Eq. © for N = 2, 
e ^ has the same form as the dual gauge field correlator 
for the case N = 1, e = 0, which exhibits a dual Higgs 
phenomenon ^lj- Thus, for N = 2, e ^ 0, — m' 2 ) 
behaves as vortices for N = 1, e = 0, i.e. it is a super- 
fluid mode arising out of superconducting condensates. A 
nonzero mAh is produced by disordering at T c \ while 
a nonzero m-A is destroyed by disordering 1 - 2 -' at T c2 . 

We have analysed the A^-color London model Eq. (J2J) 
in vortex representation Eqs. (0J and (J5J. The dual the- 
ory is given by Eqs. 10 and ©. For N — 2, we have 
performed large scale Monte Carlo simulations comput- 
ing i) critical exponents a and ii) gauge field and dual 
gauge field correlators, Hi) the corresponding masses, and 
iv) critical couplings using FSS. For + V> (2) 

we find 

one neutral low-temperature critical point at T c i, and 
one charged critical point at T c2 > T c \. For general N, a 
Higgs mass wa is generated at the highest critical tem- 
perature, followed by N — 1 neutral fixed points as the 
temperature is lowered. 

These results apply to electronic and protonic conden- 
sates in liquid metallic hydrogen under extreme pressure. 
Estimates exist for T c2 for such systems, T c2 ~ 160K 
0, and hence T c \ ~ O.IK. Hence, in addition to the 
emergence of the Meissner effect at T c2 and a corre- 
sponding divergence in the magnetic penetration length 
X ~ |1 - T/r c 2|"" /(2 " J ' A) [13, there will also be a novel 
effect, namely a low-temperature anomaly in the mag- 



as follows. Introducing \q) = (|nig 1 ' ) — m^l 2 ) and 

G^\q) = ((m« - n4 2) ) ■ (J2 2 a=1 IV^IV. ?)) we find, 
in the notation used in Eqs. ijSj and Q 
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netic penetration length A ~ 1/ttia at T c \, cf. Fig. 
due to the appearance of superfluid modes arising from 
superconducting condensates. 
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